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Abstract 

Approximate analytical chirped solitary pulse (chirped dissipative soliton) solutions of the one- 
dimensional complex cubic-quintic nonlinear Ginzburg-Landau equation are obtained. These solu- 
tions are stable and highly-accurate under condition of domination of a normal dispersion over a 
spectral dissipation. The parametric space of the solitons is three-dimensional, that makes theirs 
to be easily traceable within a whole range of the equation parameters. Scaling properties of 
the chirped dissipative solitons are highly interesting for applications in the field of high-energy 
ultrafast laser physics. 
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I. INTRODUCTION 



The complex nonlinear Ginzburg-Landau equation (CGLE) has so wide scope of applica- 
tions that the concept of "the world of the Ginzburg-Landau equation" [1| is not exaggera- 
tion. The CGLE demonstrates its effectiveness in quantum optics, modeling of Bose-Einstein 
condensation, condensate-matter physics, study of non-equilibrium phenomena, and nonlin- 
ear dynamics, quantum mechanics of self-organizing dissipative systems, and quantum field 
theory. In optics and laser physics, the CGLE provides an adequate description of mode- 
locked oscillators and pulse propagation in fibers as. 

The CGLE is multiparameter and not integrable in a general form. As a result, an 
analysis of multitude of its solutions requires extensive numerical simulations. The exact 
analytical solutions are known only for a few of cases, when they represent the solitary 
waves (dissipative solitons) and some algebraic relations on the parameters of equation are 
imposed 0, Q] ■ As a rule, one presumes some class of functional expressions to construct the 
solution. As a result of such presupposition, the solutions outside a given class are missed. 
In principle, the missed solutions can be revealed on basis of the algebraic non-perturbative 
techniques , which, nevertheless, need a lot of computer algebra. These challenges 
stimulate interest in the approximate methods of integration of the CGLE. 

The perturbative method has allowed obtaining the dissipative soliton solutions of the 
reduced and complete cubic-quintic CGLE in the limits of small and large dispersion and 
phase nonlinearity ja, B]- Another approximate approach is to reduce an infinite-dimensional 
(in terms of degrees of freedom) problem to finite-dimensional one on basis of, for instance, 
the method of moments. This allows tracing an evolution of a finite set of the trial solution 
parameters js]. 

However, some physically interesting sectors of the CGLE allow an approximation with- 
out any functional constraints imposed on the solution or/and the equation parameters. 
Moreover, dimensionality of the parametric space corresponding to such a solution can be 
reduced in comparison with the parametric space of the CGLE that makes the solution 
under consideration to be easily traceable. 

A physically important sector, which permits an approximate analysis, is represented 
by the chirped solitary pulse solutions, or the chirped dissipative solitons (CDSs) of the 
CGLE. The CDS exists in both anomalous and normal dispersion ranges ja, 10]. It is 
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121 1 and can be considered as 



Q. 



very important, that the CDS is energy scalable [8|, llO|, 111 

a model of femtosecond laser pulses with about of and over microjoule energies [13. 
Energy scalability of the CDS results from its stretching caused by a large chirp. Hence, 
the CDS with large energy has a reduced peak power that provides its stability jlol 
Simultaneously, a large chirp leads to spectral extra-broadening so that the CDS becomes 
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to be compressible down to a few of tens of femtoseconds 

The mechanism of the CDS formation is a composite balance of phase and dissipative 
effects js]. The first effect is a balance of phase contributions from the nonlinear phase 



distortion and the time-dependent phase affected by a normal dispersion [17j. That is 
possible if a soliton is chirped, but this effect alone does not provide a soliton stability. The 
CDS stability can be provided by a balance between the nonlinear gain and the spectral 
dissipation 0, 0, Q] ■ 

A large chirp of the CDS allows two main approximations: i) soliton stretching admits 
the adiabatic approximation, ii) fast phase variation allows applying the stationary phase 
method in Fourier domain [la. [la I20I] . 



dimensional CGLE Jj], |20|] and the generalized one-dimensional CGLE (i.e. the CGLE with 



As a result, the CDS of the reduced cubic-quintic one- 



an imaginary saturable law nonlinearity in terms of [19|) 16| can be represented analytically 
as the two-parametric solitary pulse solution without any restrictions on its functional form 
as well as on the equation parameters (certainly, within the scope of approximations under 
consideration, see below). 

Here the extension of this approximate technique to the complete cubic- quin tic one- 
dimensional CGLE (i.e. the CGLE with a complex parabolic law nonlinearity |l9|) is pre- 
sented. It is shown, that the CDS is the three-parametric solution with five types of the 
truncated spectral profiles: i) finger-, ii) parabolic-, and iii) flat-top, as well as iv) concave 
and v) concave-convex ones. The regions of existence and stability of the CDS are analyzed 
systematically within a whole parametric range obeying the condition of domination of a 
normal dispersion over a spectral dissipation. The obtained results are validated on the 
basis of numerical solution of the CGLE and compared with the existing results of extensive 
numerical simulations of the CGLE. 
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II. CDS OF THE CUBIC-QUINTIC CGLE 



Let the CGLE be written down in the following form 15|, |20| 



d d 2 

—a (z, t) = —era (z, t) + (a + i[3) 77-5-a (z, t)+(n— i'j) la (z, t)\ a (z, t) — (K,( + ix) \a (z, t)\ a (z, t) 
oz at 1 

(1) 

Here z is the propagation (longitudinal) coordinate, which can be the propagation distance 
in a fiber, or the cavity round-trip in a laser oscillator, for instance; t is the "transverse" 
coordinate, which can be, for instance, the local time for a propagating laser pulse 2J. The 
complex slowly-varying field amplitude a(z, t) is chosen so that \a\ 2 has a dimension of instant 
power. The first term on right-hand side of Eq. (CQ) describes an action of net-loss with the 
parameter a. In the general case, this parameter is energy- dependent (i.e. it depends on 



\a (z,t')\ dt', see 15|] ) and has to be positive to provide the vacuum stability [ i.e. the 
subcritical range of Eq. Q is under consideration]. The second term describes a spectral 
dissipation (a is the squared inverse bandwidth of spectral filter) and a dispersion [(3 is 
the dispersion coefficient]. Positivity (negativity) of (3 corresponds to normal (anomalous) 
dispersion. The third term results from a contribution of cubic nonlinearity, which is a sum 
of contributions from the nonlinear gain [or the self-amplitude modulation (SAM) defined 
by the parameter k >0] and from the self-phase modulation (SPM) with the parameter 
7. Only the focusing SPM with 7 >0 will be considered below. A higher-order (quintic) 
nonlinearity defines the fourth term in Eq. ([1]). Its real part describes the SAM saturation 
(C >0 provides stability of a desired solution against collapse), while a correction to the 
cubic SPM can be both enhancing (x >0) an d saturating (x <0). 

To find the CDS solution, let's make the traveling wave reduction of Eq.Q by means of 
the ansatz 



a [z. 



t) = ^/PVjexp[i(j)(t)-iqz}, (2) 



where P(t) is the instant power, which defines a CDS envelope; <fr(t) is the p 
the phase shift due to a slip of the carrier phase with respect to an envelope 



rase, and q is 
3]. 



Below we shall consider only the sector of CGLE, where a normal dispersion prevails over 
a spectral dissipation, that is (3 ^> a >0 [6]. This assumption is well-grounded for both 
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221 ] laser oscillators operating in the 



broadband solid-state [ll, Q, Q and fiber 
all-normal dispersion (ANDi) regime. But the numerical analysis demonstrates 16|], that 
even the case of ft > a >0 (e.g., a thin-disk solid-state oscillator [14]) can be described 
adequately in the framework of the analytical approach under consideration. 
The adiabatic approximation T ^> yfft allows obtaining from Eqs. (TIf2l) 



pn 2 = q - 7 P - X P 2 
ft dP dft\ 







+ — = K P(l-CP)-<r-aft 2 



(3) 



P(t) dt dt J 

where ft = dip (t)/dt is the instant frequency. 

Since the first equation in is quadratic in P, there are two branches of solution. 
However, it is reasonable to confine oneself to the branch, which has the limit % ^0 (this 
limit has been considered in [lsj, H])- Then, one has 



P = ~ (7 - vV + 4gx - 4 X /5ft^ 



(4) 



Since P >0 by definition, there is the maximum frequency deviation A from the carrier 
requency: A 2 = q/ ft. Thus, the second equation in (jHJ) and Eq. (jlj) lead, after some algebra 



23|, to 



dft 
~dt 



ft[4xftQ 2 -(A-^ A] 



(A — 7) A, 



(5) 



A= v/ 7 2 + 4/?x (A 2 - ft 2 ). 



The singularity points of Eq. (JHJ) impose the restrictions on the A value 23] 



A 2 



7 



(f + 1) 



2(3 + f + f)(2 + f + f ±^(c-2) 2 -16a(l + f) 

1 + f 



12 - 3c - % - 



32a 



(6) 



where three control parameters are a = cr(,l K i b = Cl/Xi an d c = aj/ftn. These three 
parameters define the parametric dimensionality of the CDS. Eqs. ( I4f6i) allow obtaining the 
CDS peak power (ft has to be equal to in Eq. (jlj) for this aim). 
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It is convenient to use the following normalizations: t' = t(n/£) y/n/aC,, A' 2 = A 2 a(/n, 
Vl' 2 = ^l 2 a(/K, P' = (P. For the dimensionless energy, one has E' = E (/c/7) yfnCJoi. 
Hereinafter, these normalization will be implied and the primes will be omitted. 

The expressions for the dimensionless quantities are 




4 (A 2 -9?) 
cb 



! ^4 (A 2 -ft 2 ) ! 



~dt 

A 2 : 



2 V V cb 

a + n 2 + ^(l-A)(l + l-A) A(A-l) 



^+A(1-A) 



16 (f + 1) 



2(f+3 + f)(f + f + 2 ±v /(c-2) 2 -16a(f + l) 

f+1 



(7) 



, 32a 
3c - 96 12 



Since the phase <p(t) of the CDS can be treated as a rapidly varying function of time in 
the limit of k <C 7, one may apply the method of stationary phase to the Fourier image of 
a(t) [20). As a result, the expression for the CDS spectral profile is 



tt (A — 1) ((A — 1) c6 + 4 (2uj 2 — A 2 )) H (A 2 — u 2 

p (u) = \e [CO) I ~ 



cA ((A -l)(c(a + b + b 2 + u 2 )+b (A 2 - to 2 )) - 2 (b + 1) (A 2 - to 2 )) ' 

(8) 

where e (to) = J dt^JP (t) exp [i<j> (t) — iut], H (x) is the Heaviside's function and one has to 
replace Q by to in A given by Eq. ([7]). 

The CDS energy can be obtained from Eq. (jSJ) by integration: E = ^p(u). This 
value can be related to the energy E* of a solution of the linearized version of Eq. ([Tl) 
through the saturated net-loss parameter a: a w 5(E/E* — 1), where 5 = da j dE\ E=E „ |l5| . 
Such a relation can be usable, for instance, to define the CDS parameters from those of a 



laser oscillator 
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Thus, the approximate technique under consideration allows representing the CDS param- 
eters, its spectral and temporal profiles as well as energy from a few of algebraic expressions 
(171IH1) . single first-order ordinary differential equation (J7J) and numerical integration of (jHJ). 
Since the CDS is three-parametric, such an approximation allows easily tracing the soliton 
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characteristics within a broad range of the CGLE parameters. It is important, that the 
absolute values of the CGLE parameters are not relevant in contrast to their relations pre- 
sented by the parameters a, b and c. This allows a unified viewpoint at the diverse systems 
obeying the CGLE 



III. MASTER DIAGRAM AND REGIONS OF THE CDS EXISTENCE 

Eq. (jTJ) demonstrates that there exist two branches of the CDS corresponding to two 
signs before square root in the expression for A. As will be shown below, such a division 
into two branches is physically meaningful. In accordance with the sign in Eq. ([7j), we shall 
denote these branches as the "positive" (+) and "negative" (— ) ones. One has note, that 
only the — branch has a limit for (, \ — > (the "Schrodinger limit") |20j. 



A. Positive branch of the CDS 

Regions of the + branch existence are shown in Fig. [T]on the plane (a-c) for the different b. 
Zero value of a corresponds to a marginally stable CDS. The existence regions are maximally 
broad in this case. The restrictions on the parameters are < c < 2, and b >0 (the SPM 
is unsaturable) or b < — c/3 — 4/3 (the SPM is saturable). One has to resemble, that the 
decrease of \b\ means a growth of contribution of the quintic SPM. The physical meaning of 
maximum c is that there is a minimum dispersion or a maximum spectral dissipation, which 
provides the CDS existence. 

A new view on the CDS results from a consideration of chirp at the soliton center (i.e. 
at t =0): 



^ = to 



b 2 

= _ a __(l_ a\ 
4 

1 + f-4 



n=o) 

(9) 



Here the chirp is normalized to k/^/3. The analysis demonstrates that the chirp becomes 
negative, when the c parameter reaches some minimum value (the lower borders of the 
hatched regions in Fig. [TJ). Such a chirp corresponds to a spike on constant background 
with lim Q < ±oo, lim P = const > 0. These solutions will be not considered hereafter 

t— *±oo t— >±oo 
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FIG. 1: Regions (hatched) of the + branch existence, b =20 (black solid curves, vertical hatch- 
ing), 0.1 (light gray dotted curves, horizontal hatching), and —2.5 (gray dashed curves, horizontal 
hatching) . 



and the chirp positivity will be admitted as the additional criterion of the CDS existence. 



This criterion agrees with the analytical results of [k| in the limit of a/f3 <C 



and K/7 <Cl. 

The appropriate zero asymptotic lim P = of the solutions analyzed in [10] exists only if 
< ~ 37/(1 + c)k < P/a. Here ip is defined as the parameter in the phase profile ansatz 



4>{t) = ^ln -y/P(t), which is used in 10]. 

As a result, there is some minimum c (i.e. maximum normal dispersion or minimum 
spectral dissipation) for the + branch (Fig. [T]). This minimum c tends to zero, when the 
positive b decreases (Fig. [Q. 

For b <0, the CDS existence range squeezes, when b approaches —2 (Fig. [[]). lib > —2, 
the positively chirped CDS has a parabolic temporal profile and lim Q = ±00, lim P = 
(r is some finite interval of local time). We will not consider such an "inverted" CDS 
hereafter. 

The a growth, if it results from the a increase, enhances the soliton stability against 
a vacuum destabilization. However, the existence regions shrink along the c parametric 
coordinate with such a growth (Fig. [1]). 
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0,00 0,03 0,06 0,09 0,12 0,15 0,18 0,21 0,24 

a 



FIG. 2: Borders of regions of the — branch existence (black solid curve for b =20, light gray solid 
curve for b =0.1, gray solid curve for b = —2.5, points for b = —2). Open circles correspond to the 



region border from 



j. The existence regions lie below the corresponding borders. 



B. Negative branch of the CDS 

Regions of the — branch existence are shown in Fig. [2] on the plane (a-c) for the different 
b. The CDS exists within the interval 0< c <2, which squeezes with a. Since this branch 
has a Schrodinger limit, such a squeezing can be obtained on the basis of the perturbative 



method 



. Then, the existence region for |6|3>lisc^2 — 4 a/ 6a/5 (open circles in Fig. 
[2]) 6|. One can see, that the approximation of [6] is quite accurate, when a«l (i.e. in the 
low-energy limit). The limiting a is denned by the hard excitation condition a ^ 1/4 (|. 
The existence region shrinks with a growing positive quintic SPM (i.e., when b >0 tends to 
zero) and stretches with a negative quintic SPM verging towards b = —2. 

There are no negative chirp solutions for this branch. There exist the positive chirp 
solutions for c >2 and b < —2, but they are the spikes on background. 

One can see that the upper (in the c direction) borders of the regions coincide for the 
positive and negative branches. This means that the branches coexist within the regions of 
their existence in the (a, b, c)-parametric space. 
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FIG. 3: The master diagram for b 3>1. There exists no CDS within the hatched region. Solid 
curve corresponds to a =0. Dashed curve divides the regions, where the + and — branches exist. 
Crosses (circles) correspond to the + (— ) branch for a =0.01. The points a and (3 correspond to 
the parameters of the numerical solutions presented in Figs. I8|9l bv open circles. The open squares 
indicate the numerically obtained stability border (k =0.047, £ =0.57). 

C. Master diagram 



Representation of the existence regions in the form of Figs. [1121 is awkward in some way, 
because the a parameter can be energy-dependent. As a result, the branches do not coexist 
as they differ in energy. It is more convenient to use a representation on the plane (E-c) for 
the different b. Such a representation will be called the master diagram. The E value can 
be easily related to the experimentally controllable parameter E* (see Section HT1). 

The master diagram for the CDS is shown in Fig. [3] for the case of vanishing quintic SPM 



(b 3>1 l24j). The solid curve shows the border of the CDS existence (a =0). Above this 



border, the vacuum of Eq. ([I]) is unstable (hatched region). The dashed curve divides the 
existence regions for the + and — branches (the branches merge along this curve). Crosses 
(circles) represent the curve along which there exists the + (— ) branch for some fixed value 
of a (so-called the isogain curve). 

The master diagram is interrelated with the existence regions in Figs. [TITS! The point of 
intersection of isogain with the dashed curve defines the maximum value of c in Figs. [T1[2l 
for the corresponding a. Since the + branch isogain has a nonzero asymptotic for E — > 00, 
there is the nonzero minimum c, which confines the + branch region for a fixed a in Fig. 
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[TJ The — branch has a zero asymptotic for E — > oo. Hence, the — branch extends down to 
c=0. 

The master diagram reveals four significant differences between the branches. The first 
one is that the — branch has lower energy than the + branch for a fixed c. The second 
difference is that the + branch isogain has nonzero asymptotic for E — > oo. In this sense, the 
+ branch is energy scalable, that is its energy growth does not require a substantial change 
of c. The — branch is not energy scalable, that is its energy growth needs a substantial 
decrease of c (e.g., owing to a dispersion growth) [if]]. The third difference is that the + 
branch verges on o =0 within a whole range of E. The fourth difference is that the — branch 
has a Schrodinger limit (, x —* 0. 

Growth of the positive quintic SPM (i.e. b 0) narrows the existence region (Fig. S]). 
This means that smaller c is required to provide the CDS existence for some E. That is, 
since the positive quintic SPM means an enhancement of the phase nonlinearity, a SPM 
enhancement has to be compensated, for instance, by a dispersion increase (c oc 1//3). One 
can see from Fig. HI that the + branch region narrows substantially with b — >0 (b >0) within 
a whole range of E. 

The situation is opposite, when the quintic SPM is negative. The existence range widens 
and a larger c (i.e. smaller dispersion) provides the CDS existence for some E. The + branch 
region widens, as well. However, it is important to remember, that the range of c, where the 
CDS with a fixed a exists, is defined by the difference between i) the point of intersection of 
the isogain with the boundary between the + and — branches and ii) the isogain asymptotic 
for E — > oo. As a result, the range of c, where some isogain exists, can be narrow in spite 
of the fact that a whole range of the + branch widens (Figs. [Tf5]) . The reversed situation, 
when a whole existence range is narrow, but the range of c for some isogain is broad, is 
possible for b >0 (Figs. ITJ3|) . 

IV. CDS PROFILE, SPECTRUM AND PARAMETERS 

Fig. [6] shows the frequency deviations and the CDS profiles relating to the + branch [see 
Eq. (0)] for the different b. One can see, that the decrease of positive b reduces a soliton 
energy (black solid vs. gray curves in Fig. ED for the fixed c and a. That agrees with Figs. 
I3ll4l where the isogain shifts towards smaller energies for a fixed c, when the positive b tends 
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FIG. 4: The master diagram for b =0.2. There exists no CDS within the hatched region. Black solid 
curve corresponds to a =0. Black dashed curve divides the regions, where the + and — branches 
exist. Gray dashed (dotted) curve corresponds to the + (— ) branch for a =0.01. The point 7 
corresponds to the parameters of the numerical solution presented in Fig. [9] by open squares. The 
open squares indicate the numerically obtained stability border (k =0.047, ( =0.57). The points 
5 and e correspond to the anal ytic al spectra presented in Fig. [12] (k =0.87 and I.57, respectively; 
(3 /a =6.25, C =0.0027, a =0.1 jl^). The point e' is the — branch counterpart of e. 
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FIG. 5: The master diagram for b =-5. There exists no CDS within the hatched region. Black solid 
curve corresponds to a =0. Black dashed curve divides the regions, where the + and — branches 
exist. Gray dashed (dotted) curve corresponds to the + (— ) branch for a =0.01. The point £ 
corresponds to the parameters of the numerical solution presented in Fig. [8] by open squares. The 
open squares indicate the numerically obtained stability border (k =0.047, (" =0.57). 
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FIG. 6: The + branch CDS profiles (solid curves) and frequency deviations (dashed curves) for 
the different b. c =1, a =0.01. 



to zero. Since a power decreases, a chirp (dfl/dt) decreases, too (black vs. gray dashed 
curves in Fig. [6]). 

In the case of b <0, the dependence of Q on t becomes "loitering" (light gray dashed 
curve in Fig. [6]). As a consequence, the CDS profile becomes flat-top (light gray solid curve 
in Fig. [6]) . The energy increases for a fixed a in agreement with a shift of the isogain towards 
larger energies in Fig. [51 

The frequency deviations and the CDS profiles for the — branch are shown in Fig. UJ Out 
of the boundary between the branches, the CDS relating to the — branch has lower energy 
and power than its + counterpart. Correspondingly, a chirp is lower, as well. Growth of the 
positive quintic SPM (i.e., b — >0) increases the soliton energy, power and chirp (black vs. 
gray curves in Fig. [TJ. 

Growth of the negative quintic SPM decreases the CDS energy, power and chirp. However, 
such a decrease is comparatively small and, therefore, it is not shown in Fig. [71 

The CDS spectra are presented in Figs. IBf9l As has been shown in Section [TT1 the spectra 
are truncated at some frequency ±A. There are the following types of spectral profiles: 
i) flat-top (solid curve in Fig. [9]), ii) convex (solid curve in Fig. [8] and dotted curve in 
Fig. [9]), hi) finger-like (dotted curve in Fig. [8]) and iv) concave (dashed curves in Figs. 
1 these typ es are widely presented in laser experiments and numerical simulations 
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One can see, that, as a rule, the CDS spectra relating to the + branch (Fig. ED are broader 
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FIG. 7: The — branch CDS profiles (solid curves) and frequency deviations (dashed curves) for 
the different b. c =1, a =0.01. 
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FIG. 8: The + branch CDS spectra for the different b. c =1, a =0.01. Solid, dotted and dashed 
curves correspond to analytical spectra. Gray circles correspond to the numerical spectrum at the 
point (3 in Fig. [3] (C = O.57, j3/a =25, k = 0.047, E =820 k^/kC, / ^/ 2 ) . Gray squares correspond 
to the numerical spectrum at the point C i n Fig- El (C = O.57, (3 /a =25, k = 0.047, E =2350 
«V></7 2 )- 



than those relating to the — branch (Fig. [9]). The cause of this difference is a smaller chirp 
for the — branch CDS. The spectrum narrows (widens) with an approach of positive b to 
zero for the + (— ) branch in accordance with a decrease (increase) of the CDS chirp. When 
the positive quintic SPM increases (b —>Q), concave spectra appear. In contrast to the model 



of 



221 ]. the source of such spectra is not the self-amplifying SAM [i.e., the negative ( in Eq. 
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a, 6 




FIG. 9: The — branch CDS spectra for the different b. c =1, a =0.01. Solid, dotted and dashed 
curves correspond to the analytical spectra. Gray circles correspond to the numerical spectra at 
the point a in Fig. [3] (two coinciding numerical profiles are defined by £ = O.I7; (3 /a =30 and 
40; k = 0.0337 and 0.0257; E =280 and E =430 ny/i^fy 2 , respectively). Gray squares 

correspond to the numerical spectrum at the point 7 in Fig. 0] (C = O.I7, (3/a =42, k = 0.0247, 
E =600 /V</7 2 ). 



(PQ)] but solely the positive quintic SPM |24|. As a result, the concave spectrum solution of 
Section 2 is stable against collapse. 

It is important to note, that a verging of b towards zero for the + branch as well as 
a transition to the — branch reduce chirp. This can violate a validity of the method of 
stationary phase (see Section [TTJ) . As a result, the spectrum edges become smooth (see 
Section El). 

As was mentioned earlier, the + branch does not vanish along the curve of o =0. This 
curve corresponds to marginal stability against a vacuum excitation and the CDS has a 
broadest spectrum here. The dependence of half-width of such a spectrum on the c parameter 
for the different b are shown in Fig. [TU1 In the absence of the quintic SPM (b 3>1), 
the dependence is symmetric relatively c =1, where the spectral width is maximum. The 
maximum A lowers (rises) and shifts towards c =0 (c =2) for the positive (negative) 6-^0 
(Fig. [T0l) . When b > —4.5, the + branch disappears for a =0 and c =1. 

Fig. [TT] demonstrates the dependencies of the spectral half- width on a for a varied b and 
a fixed c. As a result of larger energy and chirp, the + branch (solid curves) has a larger 
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FIG. 10: The CDS spectral half-widths for the + branch and the different b. a =0. 
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FIG. 11: The CDS spectral half- widths for the + (solid curves) and — (dashed curves) branches 
in dependence on a for the different b and c =1. 

spectral width, which decreases with a and the positive b verging towards zero. The region 
of the + branch existence shortens with b —+0 (also, see Fig. [1]). A negative b expands the 
+ branch region towards a larger a. However, such a region is disconnected with a =0, if 
b > —4.5 for c =1 (see the region for b = —2.5 in Fig. [1]). The existence of this minimum 
c providing the CDS with a =0 is a result of asymptotical behavior of the zero isogain in 
Fig. [5l Physically, absence of the limit a —>0 can mean that such a CDS is not able to 
develop from the vacuum of Eq. (CQ). Fig. (TT] demonstrates that the + branch disappears 
completely, when b — > —2. 
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The — branch has a lower spectral width, which increases with a (Fig. [TTJ dashed curves). 
There exists some maximum a (for a fixed c), where both branches merge. 

V. NUMERICAL SIMULATION OF THE CDS 

The above obtained approximate solution for the CDS has to be verified numerically. 
With this purpose, a symmetrized split-step Fourier method is used for numerical solving of 
Eq. (CQ). The temporal greed contains 2 16 points, and the nonlinear propagation is simulated 
in the time domain using a fourth-order Runge-Kutta method. Total propagation distance 
consists of >10 4 steps. 

The simulations demonstrate, that the necessary factor providing the CDS stability is a 
dependence of a on E. Such a dependence is chosen to be in the form presented in Section 
im (the S parameter equals to 0.5) [15j. Then, the E parameter in a master diagram can be 
easily replaced by the E* one, but the difference between E and E* is small and, therefore, 
insignificant for further consideration. 

The simulated spectra of the CDS are shown by open circles and squares in Figs. [8ll9l for 
the parameters a, b and c corresponding to the points a, [3, 7 and ( in Figs. I3"f4"f5l The 
agreement between the analytical and numerical results is perfect. Moreover, the numeri- 
cal results demonstrate that the CDS is really three-parametric and its parameters scale in 
accordance with the rules of Section [Til This means that the normalized parameters and 
profiles of the CDSs are identical for the identical sets of (a,b,c). For instance, two para- 
metric sets: i) b =20, a =0.01, /3/a =30, ( = O.I7, E* =280 k^/j 2 , k = 0.0337 [e.g., 
a 100 nJ Ti:sapphire oscillator with a =2.5 fs 2 and 7 =4.55 MW" 1 [15]]; and ii) (3/a =40, 
E* =430 Kv^C/7 2 , k = 0.0257 correspond to the single point a in Fig. [3j This is the 
— branch, and the analytical (solid curve) as well as numerical (gray open circles) profiles 
coincide in Fig. [9j A single difference between the numerical and analytical spectra is that 
the former ones have gently smoothed edges. One has note, that a scalability of the CDS 
resembles the property of a true soliton, which is a solution with not fixed parameters Q]. 

The numerically obtained stability borders are shown in Figs. 13114151 by open squares. 
The stability condition is a >0, that provides a vacuum stability. One can see, that both 
analytical and numerical borders coincide. 

It is of interest to compare the analytical results with the numerical ones presented in 
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FIG. 12: The analytical CDS spectra for the + branch at the points 5 {k =0.87; black solid curve) 
and e (k =1.57; gray dashed curve) in Fig. QJ b =0.2, j3/a =6.25, ( =0.0027, and a =0.1. 



10,ll2]- There is a difference between the parametric sectors considered in 10j,ll2| and in this 
work. The case of (3 — >0 lies beyond a validity of the analytical model under consideration, 
which requires /5>a. H (3 approaches a and then tends to zero (as well as if k prevails over 
7), the spectrum edges become smooth [3,[l6] rather than truncated. The scaling rules of 



Section [IT] and the requirement of c <2 can get broken in this case 25]. 

Nevertheless, i) strong scalability of E with j3, as well as both ii) existence of maximum 
and iii) minimum k providing a stable soliton suggest that the solutions analyzed in fiol \v* 
belong to the + branch (here we consider only normal dispersions). In conformity with |l0| . 
the stable CDS exists within the region of normal dispersion ((3 >0) for both positive and 
negative x (Figs- HE]). A fast disappearance of the CDS with the increase of b < [h]] is 
the characteristic feature of the + branch (Fig. [I]). 

As expected, the CDS of ljj belongs to the + branch (the points 5 and e in Fig. HI). 
The corresponding analytical spectra are shown in Fig. [12] [p{uj) for the black solid curve 
is re-scaled for convenience]. Both analytical spectra match with the numerical ones in Fig. 
3 of [12] with the exception of the smoothed edges for the latter owing to k > 7. Such a 
smoothing enhances for the — branch (the point e' in Fig. H]) because a chirp is lower for 
this branch. 

When k < 7 , the spectrum is concave (black solid curve in Fig. fT2l) like that in Fig. [8] 
for b =0.2. When k exceeds 7, the new type of a spectral shape appears: the concave-convex 
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one (gray dashed curve in Fig. [T2|) . Such spectra have been studied numerically in 12j for 
k > 7. 



VI. CONCLUSION 



In conclusion, approximate chirped solitary pulse solutions of the cubic-quintic nonlinear 
CGLE have been constructed analytically under condition of domination of a dispersion 
over a spectral dissipation. The solutions are three-parametric and easily traceable within a 
whole parametric space, which has been represented in the form of the master diagrams. The 
solutions are divided into two branches, which differ in their energies and scaling properties. 
It is found, that the chirped dissipative solitons under consideration have truncated spectra 
with the concave, convex and concave-convex tops. Numerical analysis and comparisons with 
the existing results have demonstrated, that the approximate analytical solutions are stable 
and highly- accurate. The obtained results are of interest, in particular, for a development of 
both solid-state and fiber laser oscillators aimed to a generation of femtosecond pulses with 
over- micro joule energy. 
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